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We obtain the phase diagram of spin-imbalanced interacting Fermi gases from measurements of 
density profiles of (> Li atoms in a harmonic trap. These results agree with, and extend, previous 
experimental measurements. Measurements of the critical polarization at which the balanced su¬ 
perfluid core vanishes generally agree with previous experimental results and with quantum Monte 
Carlo (QMC) calculations in the BCS and unitary regimes. We disagree with the QMC results in 
the BEC regime, however, where the measured critical polarizations are greater than theoretically 
predicted. We also measure the equation of state in the crossover regime for a gas with equal 
numbers of the two fermion spin states. 
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Strongly interacting Fermi gases are found in a 
variety of settings, including superfluid 3 He, quark 
matter, superconducting materials, and ultracold 
atomic gases CH3j. The properties of such sys¬ 
tems, including the nature of any superfluid or su¬ 
perconducting order, strongly depend on the inter¬ 
actions between particles. At sufficiently low tem¬ 
peratures the short-range interaction between oppo¬ 
site spin atomic fermions may be characterized by 
the parameter 1 /kpa, where kp is the Fermi mo¬ 
mentum, and a is the s-wave scattering length. For 
1 /kpa > 1, opposite spins may form tightly bound 
bosonic pairs which repel each other, thus creating a 
Bose-Einstein condensate (BEC) of molecules. For 
weaker attraction, where 1/kpa < —1, an ultra¬ 
cold atomic gas may form a conventional Bardeen- 
Cooper-Schrieffer (BCS) superfluid of loosely bound 
pairs. In between these extremes is the unitarity 
regime, —1 < 1/kpa < 1, corresponding to resonant 
two-body interactions. This BEC-BCS crossover 
has been studied extensively over the past decade 
in the context of ultracold atomic Fermi gases m- 
When the two spin states have equal populations, 
the crossover between the BEC and BCS limits has 
no phase transitions as a function of l/kpa. Addi¬ 
tional phases can appear, however, when an effective 
magnetic field couples to the spin-!/2 fermions, favor¬ 
ing an imbalance (or polarization) in the number of 
fermions in each spin state 00 - In thin-film elec¬ 
tronic superconductors, such a coupling can come 
from a real in-plane magnetic field [9] . In the present 
setting of cold atomic gases, this imbalance is accom¬ 
plished by creating unequal populations of the two 
hyperfine levels comprising the pseudo-spin-*/2 sys¬ 
tem. In the BCS regime, a sufficiently large chemical 
potential difference, known as the Chandrasekhar- 
Clogston (CC) limit [101 ITT], will suppress pairing. 
A spin imbalance can be accommodated in the BEC 


regime, however, resulting in a Bose-Fermi mixture 
that remains a superfluid. The exotic Fulde-Ferrell- 
Larkin-Ovchinnikov (FFLO) state, featuring pairs 
with non-zero momentum, has been proposed as the 
ground state of a spin-imbalanced superconductor 
under certain conditions Qjug. There have been 
no definitive observations of FFLO superconductiv¬ 
ity, but an experiment on spin-imbalanced fermions 
confined to one dimension has produced a phase di¬ 
agram with a large polarized region consistent with 
FFLO M- 

Neglecting any exotic superfluid phases (such as 
the FFLO), the phase diagram of the three dimen¬ 
sional (3D) spin-imbalanced Fermi gas as a func¬ 
tion of interaction strength and polarization, ex¬ 
hibits four phases [15]: (i) fully polarized, non¬ 

interacting normal (Npp), (ii) partially polarized 
normal (Npp), (iii) partially polarized superfluid 
(SFp), and (iv) unpolarized superfluid (SF 0 ) T6J. 
Additionally, constraining the system to fixed parti¬ 
cle number leads to regions of phase-separated mix¬ 
tures of these phases. The local polarization, defined 
as the effective magnetization divided by the density, 
p = (n-f — n^)/(n^ + n 4 .), vanishes in the SF 0 phase, 
p = 1 in the Npp phase, and 0 < p < 1 in the Npp 
and SFp phases. The majority and minority species 
are defined by n-j- > rip 

Experimentally, atoms are generally trapped in 
potentials resulting in inhomogeneous density distri¬ 
butions. In the local density approximation (LDA) 
the local state of the gas is determined by its local 
chemical potential, so the density profiles can reveal 
transitions between phases. Observations of phase 
separation in spin-imbalanced Fermi gases were ob¬ 
tained at Rice and MIT by direct in-situ imaging of 
the density distributions QMS], and by imaging the 
distributions in time-of-flight m- The distributions 
in the Rice experiment were out-of-equilibrium due 


2 


to an evaporative depolarization mechanism at work 
in their highly elongated confining potential nn m, 
and could not, therefore, be compared with distribu¬ 
tions calculated assuming equilibrium. Density pro¬ 
files obtained by the MIT group at unitarity and on 
the BEC side of resonance [ 23 ] agree quantitatively 
with the theory of Bertaina and Giorgini (BG) com¬ 
puted using quantum Monte Carlo (QMC) and the 
LDA [ 23 ]. In the unitary regime, these profiles con¬ 
tain a jump in the local polarization p that indicates 
a first-order phase transition between the superfluid 
and the normal phases [ 23 ] [ 23 ]. The ENS group mea¬ 
sured thermodynamic properties of the imbalanced 
gas by extracting the equation of state from doubly- 
integrated density profiles [25] • In this paper, we 
report new measurements of the density profiles of 
Fermi gases for —1 < 1 /fcj?a< 2 and use these mea¬ 
surements to better constrain the low-temperature 
phase diagram. These measurements largely confirm 
the results of previous investigations and extend the 
range of interactions studied. 

Our method for producing an imbalanced degener¬ 
ate gas in the lowest two hyperfine states of 6 Li, F = 
1/2, m F = 1/2 (|t)) and F = 1/2, m F = -1/2 (|t», 
has been discussed previously in detail mmm- 
In brief, we sympathetically cool 6 Li with 7 Li in an 
Ioffe-Pritchard magnetic trap, then load the 6 Li into 
a single-beam optical dipole trap formed by a fo¬ 
cused infrared laser beam. We control the spin im¬ 
balance by varying the power of an adiabatic RF 
transfer from |f) to |t) at a field of 835 G. After the 
RF transfer, we evaporatively cool the cloud in the 
single-beam trap by reducing its depth. We evapo¬ 
rate at 835 G to study interactions on the BCS side 
of the broad Feshbach resonance at 832 G HU ]23, 
while for fields on the BEC side of resonance we 
quickly ramp the field to 765 G before evaporation. 
After evaporation, atoms are loaded into the final 
trap formed by two focused infrared laser beams 
crossing at right angles while the single-beam trap is 
slowly (100 ms) ramped off. The crossed beams each 
have 1/e 2 radii of 55 (im x 235 pin. resulting in an 
ellipsoidal crossed-beam trap with a measured axial 
frequency of u> z / 2-7T = 78 Hz and measured radial fre¬ 
quencies of lo x /2it = 248 Hz and u y /2n = 274 Hz, at 
a trap depth of 1.5 pK. The number of |f) atoms, A^, 
is typically around 2 x 10 5 , and varies by about 10% 
shot-to-shot. The cloud polarization P = 
varies from shot-to-shot by about 30% for a given 
RF power, so data must be post-selected using the 
measured P. After loading into the crossed-beam 
trap, we ramp the magnetic field to its final value B 
at a rate from 0.4-2.0 G/ms; the final bias field has 
an uncertainty of 2 G. 

We use in situ phase-contrast polarization imag¬ 
ing, described previously [14, [28] to record the spa¬ 


tial distribution of the trapped atoms. The probe 
beam propagates perpendicular to the bias magnetic 
field, which is parallel to the axial trap direction. 
The column densities, n^^x', z), of each spin state 
are extracted from two images taken within sev¬ 
eral ps of each other at different probe detunings. 
Here the imaging plane ( xz ) is rotated 30° from 
the (x, z) plane defined by the trap potential. To 
improve the signal-to-noise ratio, we fit n c to find 
the cloud center (x r = 0,z = 0), then average the 
four quadrants to obtain the column density distri¬ 
butions of the majority, n c( p{x',z) = ^[n c ^(x', z) + 
n c ^(-x', z) + n c ^(x', -z) + n c ^(-x', —z)], minority 
n cq i(x', z), and their difference, which is related to 
the spin density. The top row of Fig. |T] shows the 
average of these column densities for several exper¬ 
imental realizations with fixed parameters, for sev¬ 
eral values of B and P. The majority and minor¬ 
ity cloud radii, R^ and R j,, are obtained from axial 
cuts of the column densities for each experimental 
run, then averaged over several runs. We also deter¬ 
mine the radius, R s , where the spin column density, 
n cq ^ — n cg 4 .) is maximum (the ‘cusp’)- Within the 
LDA, a cusp with a discontinuous derivative would 
indicate the location of a first-order phase transition 
for a uniform gas. These mean radii are indicated 
by the vertical lines in Fig. [T] 

We reconstruct the density distributions n-pj.^) 
using inverse Abel transforms of the averaged n cq -fp. 
The bottom row in Fig. [I] shows axial cuts of these 
density distributions. The SF 0 core radius, R C1 is 
the radius at which the spin density first rises above 
zero. The mean radii for several experimental re¬ 
alizations are indicated by the vertical lines in the 
bottom row of Fig. [T] Experimentally, we determine 
R c by finding where the spin density first rises above 
the background spin density noise, which is the stan¬ 
dard deviation of the spin density for z > R^. To re¬ 
duce bias toward obtaining smaller values of R c due 
to noise, we smooth the profiles with a 7-pixel-wide 
Hann window before computing R c . We also confirm 
our determination of R c by fitting the spin density 
profiles near R c with a function that increases lin¬ 
early from 0 for z > R c ‘, the fit results are consistent 
to within shot-to-shot variation. 

Temperatures are measured by fitting the ferro¬ 
magnetic wings of n c | for clouds with high P to 
non-interacting Thomas-Fermi distributions. We 
find that for B > 743 G the fitted temperature 
T < 0.08Tp, where T F ss 1.5 pK is the Fermi tem¬ 
perature of N-f non-interacting atoms. For lower val¬ 
ues of B , however, we measure higher temperatures, 
which are likely a result of heating from inelastic 
molecular decay collisions. At B = 725 G, for ex¬ 
ample, we find T « 0.11 Tp. 

The boundary locations -R-j-, R±, R Sl and R c , are 
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B = 736 G.P = 0.60 B = 780 G, P = 0.73 B = 832 G, P = 0.50 B = 832 G, P = 0.90 B = 890 G, P = 0.40 



0 20 40 60 80 100 
Radius (pm) 


FIG. 1. (Color online) (a)-(e) Axial cuts of the quadrant-averaged column density n cq -|-,.|,(z) and (f)-(j) their corre¬ 
sponding density profiles Values of B and P are indicated for each column, and the corresponding values 

of 1/fc.Fta are: (a, f) 1.6; (b, g) 0.6; (c, h), (d, i) 0; (e, j) —0.4. The uncertainty in l/fcF-|-a is as large as 0.07 
based on a combination of a 10% systematic uncertainty and shot-to-shot variation in N-f, 3% uncertainty in the trap 
frequencies, and 2 G uncertainty in the bias magnetic field. Each plot is an average of 3-9 experimental realizations 
that have P within a range AP = 0.02 centered on the given value. Black, blue and red curves correspond to the 
majority spin |'|'), minority spin |4_), and their difference, respectively. In the upper row, the blue and black vertical 
lines indicate the mean of the minority and majority edges, Ri and R-f, respectively, and the purple vertical lines 
indicate the radius of maximum column density difference R s . In the lower row, the vertical green lines indicate 
the mean boundary of the SFo core, R c . For each vertical line, the standard error of the mean is indicated by the 
line’s thickness. We estimate a systematic uncertainty in the radii of 4 fim, dominated by the resolution limit of our 
imaging system. For (d, i), P > P c so that R c = 0, and R s is not meaningful. 


plotted as functions of P in Fig. [2] for several dif¬ 
ferent interaction strengths, l//cp-|-a, ranging from 
the BEC to the BCS regimes. These boundary radii 
are normalized by R z = , 

the axial Thomas-Fermi radius for a non-interacting 
gas with iVf atoms, where a z = ( h/mu > z ) 1 / 2 is the 
axial harmonic oscillator length. The interaction 
strength is determined from fcpf = (48A r |) 1 / 6 /aho 
and a = a(B) [27], where dho = (fi 3 /m 3 WjW I w s ) 1 / 6 
is the mean harmonic oscillator length, and B is the 
bias magnetic field. For a given P, the systematic 
variation in with P produces up to a factor of 
1.2 variation in 1 /kp^u. Due to this variation, ex¬ 
periments at a given field trace out the P-A/k,F\a 
phase diagram along non-vertical lines. To account 
for day-to-day variation in trap frequencies, we scale 
R z for all the data at a given B so that R^/R z goes 
to 1 as P goes to 1 -this variation is less than 5%. 

The radii plotted in Fig. [2] provide detailed in¬ 
formation about the phases of trapped imbalanced 


Fermi gases as a function of the imposed popula¬ 
tion imbalance. One common feature is the exis¬ 
tence of a balanced SFo core with radius R c that 
decreases with increasing P until it vanishes at a 
critical cloud polarization, P c . To extract P c we fit 
R C (P) for each field, shown by the green data points 
in Fig. [2] to an empirical function which vanishes 
as (P c — P) 1 / 2 for P < P c . The results are indi¬ 
cated by vertical green arrows in Fig. [2] At unitar- 
ity, we measure P c = 0.79(4), where the error bar 
accounts for the uncertainty in measuring P for a 
single cloud as well as systematic uncertainty in the 
best fit parameters. This result is in good agreement 
with previous measurements giving P c = 0.77 [18] , 
0.76(3) [ 22 ]: and 0.75 [ 251 , and with theoretical pre¬ 
dictions of P c = 0.77 [MUSI, all slightly higher than 
an initial measurement of P c = 0.70(3) [ 201 . 

The work of BG, following earlier calculations of 
Pilati and Giorgini m, involved calculating the 
phase diagram of trapped Fermi gases by combining 
the LDA with fits to QMC calculations to character- 
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B = 725 G, 1//fpta = 2.0-1.6 B = 743 G, 1//t Ft a = 1.4-1.2 B = 765 G, 1//( Ft a = 0.88-0 80 B = 780 G, l/fr^a = 0.62 


B = 800 G,1//t Ft a = 0.34 
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FIG. 2. (Color online) (a)-(e), (k)-(o) Radii extracted from density and column density profiles at several interaction 
strengths: the majority Rf, (A), minority Rf (Y), cusp R a (■), and SFo core R c (•) radii as functions of P, scaled 
by the axial Thomas-Fermi radius, R z , of a non-interacting Fermi gas with Nf particles, (f)—(j), (p)-(t) Local 
polarization at the cloud center, po (•), and at R a , p a (■). Each data point is the average of several realizations of 
the experiment, binned with width AP = 0.02. Some of the phase boundaries R a and R c could not be identified for 
small P due to poor signal-to-noise, and for high P there is no identifiable R a . In these instances, the data points 
are omitted. The values for 1/kFfa have uncertainty less than 0.07, resulting from 10% systematic uncertainty and 
shot-to-shot variation in Nf, 3% uncertainty in the trap frequencies, and 2 G uncertainty in the bias magnetic field. 
However, due to systematic variation of Nf with P, 1/kFfa varies with P for a given B, particularly in the deep BCS 
and BEC regimes. In these cases, we list a range of values from l/kpfO, at P = 0 to the value at P = 1, otherwise, we 
list the mean value. For each interaction strength, R c decreases as P increases, until vanishing at P c (green arrow), 
which we determined with a fit (see text). We also fit po to determine P c (orange arrows), as described in the text. 


ize the ground-state energies of the strongly interact¬ 
ing balanced SFo phase and the partially polarized 
normal phase Npp. The ground-state energy of the 
SF p was taken to consist of contributions from a bal¬ 
anced superfluid of pairs (given by the SFo equation 
of state), a noninteracting Fermi gas of the excess 
|t) spins, and a leading-order interaction between 
|t) spins and pairs characterized by the atom-pair 
scattering length at/ = 1.18a. This characterization 
of the SFp state was found by Pilati and Giorgini 
to agree quite well with their QMC calculations. In 
addition, we have repeated the BG calculations in¬ 
cluding additional terms in the expression for the 
ground-state energy. The theory of BG includes 
an interaction between Cooper pairs, with density 
oc n±, and excess |f) spins, with density oc nf — nq, 
resulting in an interaction strength oc n±{rif — 74). 


Work by Alzetto and Leyronas has found a higher- 
order correction with strength oc rif{nf— rzq) 4 / 3 [311 . 
However, we find that including this correction (con¬ 
tained in Eq. (53) of Ref. [3T|) within the BG for¬ 
malism does not appreciably alter the value of P c for 
any interaction strength. Thus, we expect the BG 
result for P c as a function of interaction strength to 
be a robust theoretical prediction that we can test 
with our measurements. 

We also determine P c by finding the value of P 
where the local polarization at the center of the 
cloud, po = p{z = 0), Erst rises above zero by 
fitting po(P) to a function that increases with P 
for P > P c . For 1/kFfa > 0.5, where we find a 
continuous SFo-SFp phase boundary, we assumed 
that Po(P) increases with a sum of terms going like 
(P —P c ) 3 / 2 and (P —P c ) 5 / 2 . This form is motivated 
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by the mean-field result for the magnetization M 
vs. chemical potential difference in the SFp state of 
a 3D Fermi gas [8], 

2 m 3 / 2 / ,_ \ 3/2 

0{h ~ hA 

(!) 

with m the atom mass, h the chemical potential dif¬ 
ference, fj, the chemical potential and A the local 
pairing amplitude. As seen by the presence of the 
Heaviside step function 0(h — h c ), the magnetiza¬ 
tion is nonzero only for h > h c = -\/|/i| 2 + A 2 and, 
close to h c the onset of M is a sum of terms going as 
( h—h c ) 3 / 2 and (h — h c ) 5 ^ 2 . If we furthermore assume 
that, at low P , the cloud polarization scales linearly 
with h, then we have justified our assumed form for 
the behavior of po(P), allowing us to extract P c . 

Away from the deep BEC regime, for l/kp^a < 
0.5, where we find a first-order phase transition 
SFo—>Npp [51], we fit po(P) with a function that 
is linear in (P — P c ), the expected magnetization 
for a Pauli paramagnetic phase. The values of P c 
obtained from these fits are indicated by vertical or¬ 
ange arrows in Fig. [2j While our two methods should 
ideally produce the same P c , they differ slightly be¬ 
cause we only consider non-negative radii, which 
leads to slight overestimates of R c near P c when av¬ 
eraging several profiles. The magnitude of this effect 
is smaller than the uncertainty in determining P. 
Furthermore, due to noise in the density profiles, we 
cannot distinguish an SFo core from an SFp phase 
with p < 0.03. 

The dependence of the critical polarization P c on 
1 /kp-\a determined by both methods is shown in 
Fig. p| P c reaches a maximum near \/kp\a = 0.7 
and decreases as the interactions are tuned in either 
direction. Our measured values of P c agree with 
the values from the MIT [T5] [20] and ENS 1251129| 
groups for 1/fcp^a < 0.75, where our measurement 
ranges overlap. Our measurements also agree with 
the zero-temperature BG theory 1241 in this regime. 

According to theory, P c begins to drop for 
l/fep-f-a > 0.7, as the BCS pairs transition to more 
tightly bound, bosonic molecules m- As 1/kFfO, 
increases, the superfluid becomes more bosonic in 
character than fermionic, and since the bosonic su¬ 
perfluid can accommodate free fermions, the SFo 
core begins to vanish. For \/kp\a < 0.7, the transi¬ 
tion from SF 0 to Npp is predicted to be first order, 
while for l/kp^a > 0.7 the transition from SFo to 
SFp is continuous m ■ In this Bose-Fermi regime, 
we observe critical polarizations for loss of the unpo¬ 
larized core to be somewhat higher than predicted 
by BG [23]. It is unlikely that this discrepancy is due 
to the elevated temperatures we obtain in the BEC 
regime, since P c is expected to decrease with increas- 



FIG. 3. (Color online) Critical polarization of a trapped 
gas, P c , as a function of the interaction parameter 
l/kF-[CL in the BEC-BCS crossover. An unpolarized su¬ 
perfluid core exists for P < P c . The green points are the 
value of P at which the SFo core radius vanishes based on 
the fits described in Fig. [2] The orange points show the 
value of P above which the polarization at the center of 
the cloud is nonzero based on fitting to the appropriate 
function (see text). Vertical error bars include the uncer¬ 
tainty in determining P of 0.03, measured by preparing a 
series of known balanced clouds and finding the variation 
of P, as well as uncertainties in fitted parameters. From 
unitarity to the BCS side, our results agree with previ¬ 
ous experimental results from MIT [181 . i 20] , (open and 
closed red circles, respectively) and from ENS (black 
line) [25], as well as with the theory of BG (red dashed 
line) [ 23 ]. For 1/kFfO, > 0.7, we find P c to be higher 
than predicted by BG. 


ing T [32]. Thus, the effect of finite T is to diminish 
the SF 0 phase in favor of the SFp phase, while we 
actually observe a more robust SFo phase. Another 
possible explanation is that the discrepancy arises 
from the experimental challenge of observing a small 
central polarization that increases from zero contin¬ 
uously with increasing P, rather than as a first-order 
jump, as in the BCS regime. 

As we have discussed, the critical polarization P c 
indicates where the balanced superfluid core of a 
trapped gas disappears. The data also reveal infor¬ 
mation about the uniform density phase diagram, 
assuming the LDA holds. To study these phase 
boundaries, we measure the local polarization p s at 
the radius of maximum column density difference, 
p s = p(R s ), for each cloud. According to the LDA, 
jumps in the atom density as a function of chemical 
potential in a uniform imbalanced gas lead to jumps 
in the density profile in the trapped gas. These 
jumps occur at radii of maximum column density 
difference, implying that p s can indicate the critical 
polarization for a first-order phase transition in the 
uniform system [241 . 

In Fig. [2] the second and fourth rows show the de- 
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FIG. 4. (Color online) Critical local polarization of a 
homogeneous imbalanced gas, extracted by finding the 
polarization p c at the cusp location (where a first order 
phase transition occurs), as a function of the interaction 
parameter l/fept® in the BEC-BCS crossover. Vertical 
error bars reflect the standard deviation of p a (P) for P 
within AP = ±0.05 of P c . Our results agree with theory, 
shown as a green line m, though we find somewhat 
higher p c than previous experimental results, indicated 
by red points [23] . 


pendence of p s on the cloud polarization P. We can 
identify three distinct regimes showing qualitatively 
different behavior. First, in the deep BEC regime, 
1/fcj^a > 0.7, we observe that p s increases to 1 as 
P goes to 1. This behavior indicates that, in this 
regime, p s does not measure the position of a uni¬ 
form system phase boundary within the LDA, but 
is instead simply a local maximum of p within an 
SFp phase. In this coupling range, therefore, the 
critical polarization for the superfluid transition of 
a uniform gas is p c = 1 [31] ■ 

Near the unitary region, for l/hp^a < 0.3, p s is 
seen to be approximately independent of P for a 
wide range of P (see, e.g., Fig. |2 Jd) . The presence of 
the plateau indicates that the LDA holds, and that 
the point of maximum column density difference in¬ 
deed represents a jump in p and a corresponding 
phase transition in the uniform system at p c between 
SF 0 and Npp phases. We take p c to be the mean 
value of the plateau for P < P c . 

Finally, in between these two regimes, for 0.3 < 
l/fcp^a < 0.7, we find that p s increases monotoni- 
cally with P, but that p s = 1 is never reached. This 
is the regime, predicted by BG, in which there is an 
SFp phase, but only for sufficiently small P. Here, 
we take p c to be the asymptotic value of p s evaluated 
at P c . The values of p c extracted in these regimes 
are plotted in the phase diagram (Fig. |4|, and show 
excellent agreement with QMC calculations m- 
Possible evidence for finite temperature, and per¬ 
haps finite imaging resolution, is the absence of clear 


jumps in the minority density profiles shown in Fig. [T| 
for the unitarity and BCS regimes where a first order 
transition between SFo and Npp phases is expected. 
Systematic effects are also evident in the phase dia¬ 
grams of Fig. [2j In the unitarity/BCS regimes, R s 
should correspond to R c , whereas in the BEC regime 
for 1/fcp^a > 1, R„ should correspond to R±, since 
the transition is between SFp and Npp phases [23] , 
While the predicted trends are observable in the data 
the agreement is not exact. 

Finally, the Equation of State (EoS) of a bal¬ 
anced (P = 0) gas is given by £(1 /kp^a) = (Psf — 
^-Eb)/{^NEp), where E$f is the ground state en¬ 
ergy of the superfluid, Ep is the Fermi energy, N is 
the total number of atoms, and Eb = — h 2 /ma 2 is 
the binding energy for a molecular pair when a > 
0 US US] • For a harmonically trapped gas at unitar¬ 
ity, the EoS can be rewritten as £(0) = ( R p /R z ) 4 , 
where R p is the radius of the superfluid core and 
R z is the Thomas-Fermi radius of a non-interacting 
gas with N-f particles [331131]. Although (R p /R z ) 4 
only approximates the EoS away from unitarity, we 
nonetheless present our measurements of this quan¬ 
tity in Fig. [3] and compare them with previous mea¬ 
surements and theoretical calculations of the EoS. In 



FIG. 5. (Color online) The Equation of State (EoS) 
£(l/&Ft a ) = (E S f — ^Eb)/{^NEF) for an unpolarized 
gas. The data points show (R p /R z ) 4 , where R p is the 
superfluid core radius, and R, is the axial Thomas-Fermi 
radius of a non-interacting Fermi gas with N-f atoms. At 
unitarity, and on the BCS side of resonance, we take 
R p — Ri at P = 0, while on the BEC side, we deter¬ 
mine R p by extrapolating R s to P = 0. Since N )- varies 
with P for a given field, the values of 1/kF^a for P = 0 
differ slightly from those at P = P c , as in Figs. [3] and [4] 
At unitarity, we find £(0) = 0.39(3) in agreement with 
a previous measurement from MIT [4Tj (red point). Al¬ 
though ( Rp/R z ) 4 only approximates the EoS away from 
unitarity, our results agree with theoretical predictions of 
the EoS [16] (dashed green line), and with experimental 
results from ENS [25] (black line). 
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the BEC regime, we fit the column density profiles 
to a sum of Thomas-Fermi and Gaussian distribu¬ 
tions. Since we find that the Thomas-Fermi radius 
corresponds to R s for low P, we find R p by linearly 
extrapolating R s to P = 0. At unitarity and in the 
BCS regime, where the superfluid is unpolarized, we 
take R p = R± for data with P = 0. At unitarity, we 
find £(0) = 0.39(3), in good agreement with theo¬ 
retical calculations of the Bertsch parameter [351110] 
and recent measurements Esnn]. 

In conclusion, we have measured density profiles 
of spin-inrbalanced Fermi gases across the BEC-BCS 
crossover. From these profiles, we determined the 
critical polarization for both harmonically trapped 
and uniform gases above which the balanced super- 
fluid phase SFo is suppressed. The agreement with 
previous measurements and QMC theory is generally 
good, although we find a more robust SFo core in the 
BEC regime than predicted by theory. Although this 
discrepancy may be explained by very small polar¬ 


izations that are difficult to detect, the data show 
that we are able to resolve po as small as 0.03. It 
may also be possible that small adjustments to the 
theory could result in relatively large changes to P c . 
Finally, we have measured the equation of state in 
the crossover regime, which is consistent with theory 
from the BCS to the BEC regimes. 
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